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Abstract 

Quarks whose left- and right-handed chiral components are both singlets 
with respect to the SU(2) weak-isospin gauge group, offer interesting physics 
possibilities beyond the Standard Model (SM) already studied in many con- 
texts. We here address some further aspects. We first collect and update 
the constraints from present data on their masses and mixings with conven- 
tional quarks. We discuss possible effects on b — ► s~f and Z — > bb decays 
and give fresh illustrations of CP asymmetries in B° decays differing dramat- 
ically from SM expectations. We analyse singlet effects in grand unification 
scenarios: d-type singlets are most economically introduced in 5 + 5* mul- 
tiplets of SU(5), with up to three generations, preserving gauge coupling 
unification with perturbative values up to the GUT scale; u-type singlets 
can arise in 10 + 10* multiplets of SU(5) with at most one light generation. 
With extra matter multiplets the gauge couplings are bigger; we give the 
two-loop evolution equations including exotic multiplets and a possible extra 



U(l) symmetry. Two-loop effects can become important, threatening unifi- 
cation (modulo threshold effects), perturbativity and asymptotic freedom of 
03. In the Yukawa sector, top-quark fixed-point behaviour is preserved and 
singlet-quark couplings have infrared fixed points too, but unification of b and 
r couplings is not possible in a three-generation Eq model. 
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I. INTRODUCTION 



In addition to the established three generations of quarks in the Standard Model (SM), 
the possible existence of exotic singlet quarks (whose left and right chiral components are 
both singlets with respect to the SU(2) weak isospin gauge group) has been raised in various 
contexts. It was once questioned whether the b quark might be such a singlet, with no 
doublet partner t One charge — | singlet quark appears naturally in each 27-plet fermion 
generation of E 6 Grand Unification Theories (GUTs) Charge | singlet quarks have 

been variously motivated, as part of a new mass mechanism for top quarks [H or as part of 
a new supersymmetric gauge model with natural baryon-number conservation ||. If they 
exist, both kinds of singlet quarks can be produced via their strong and electroweak gauge 
couplings; mixing with standard quarks then allows the mixed mass eigenstates to decay 
via charged currents (CC) or neutral currents (NC) to lighter quarks q plus W or Z |§,§j,[7|, 
and also via Yukawa couplings to q plus Higgs bosons H P|§ . Singlet quark production and 
decay can therefore give characteristic new signals and modifications of old signals, discussed 
in the literature |]3|-p|,|7|,|8|JT0|-[16| . Possible indirect consequences of singlet-quark mixing for 
flavor- changing neutral currents (FCNC), flavor-diagonal neutral currents (FDNC) and CP 
violation have also been considered P,[TT|-p6[| 

In the present paper we address some further aspects of singlet quark physics. We first 
collect and update the direct and indirect constraints on masses and singlet-doublet mixing 
from present data, illustrating possible effects on b — > 37 and Z — > bb decays and on CP 
asymmetries in neutral B decays. We then analyse the impact of Q = — | and Q — | singlet 
quarks on the renormalization group equations (RGE), on the unification and perturbativity 
of gauge and Yukawa couplings, and on the exotic matter multiplets in GUT scenarios. 

Section II introduces our notation and lists general basic properties of singlet quark 
couplings and mixings with SM quarks. Section III addresses the 4x4 mixing matrix, 
arising when one singlet mixes with three SM quarks, and extracts the full set of unitar- 
ity constraints based on present limits on the CKM submatrix. Section IV discusses the 
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constraints implied by the absence of identifiable signals from singlet-quark production and 
decay at present e + e~ and pp colliders. Section V considers tree and box diagram contribu- 
tions to neutral meson-antimeson oscillations and the indirect constraints on singlet-quark 
mixing from present data. Section VI addresses indirect constraints from FCNC and FDNC 
decays, including a new more stringent measurement of K L — > and weak bounds from 

B°,D° — > limits; the topical cases b — > 57 and Z — > bb are discussed here. The global 

FDNC constraints are comprehensive enough to have useful repercussions via unitarity, for 
d-type singlet mixing. Section VII discusses CP asymmetries in neutral B decays, with 
new illustrations of how rf-type singlet mixing can give dramatic changes from SM expec- 
tations. Section VIII, our major new contribution, analyses the possible roles of singlet 
quarks in GUT scenarios. We show that rf-type singlets are most economically introduced 
in 5 + 5* multiplets of SU(5), with up to three generations, preserving gauge coupling 
unification and perturbativity up to the GUT scale; w-type singlets can arise in 10 + 10* 
multiplets of SU (5) with at most one light generation. The presence of extra matter mul- 
tiplets makes the gauge couplings bigger and two-loop effects potentially more important. 
We give the two-loop evolution equations, including the effects of exotic matter multiplets 
and a possible additional U(l)' gauge coupling, and show that two-loop effects can threaten 
not only unification (where threshold effects may partly compensate) but also perturbativity 
and asymptotic freedom of a 3 at large scales. In the Yukawa sector, top-quark fixed-point 
behaviour is preserved and singlet-quark couplings have infrared fixed points too, but unifi- 
cation of b and r couplings is not possible in a three-generation Eq model. Finally, Section 
IX summarizes our conclusions while Appendices A and B contain some technical details. 

II. BASIC PROPERTIES AND NOTATION 

We shall generally denote singlet quarks by the symbol x, and SM quarks by q. More 
specifically, Xd denotes a generic charge — | singlet and x u implies charge |. The weak 
isospin T 3 and hypercharge \Y of the left and right chiral components, characterizing their 



4 



SU(2) and U(l) gauge couplings, contrast with SM assignments as follows: 
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where Q = T 3 + |Y is the electric charge. The vector and axial couplings to Z are 

u d x u x d 

9v \ — fsin 2 ^ — i + lsin 2 ^ — | sin 2 ^ ^sin 2 ^ 

1 1 
9a —j 4 

where 9y/ is the Weinberg angle. Both SM and singlet quarks are color triplets and have 
the same couplings to gluons g. Hence singlet quarks have pure vector gauge couplings to 
g>,7, Z (and zero coupling to W); they are sometimes called "vector-like" or more precisely 
"vector-singlet" quarks. They do not contribute to chiral anomalies. 

(a) 2 x 2 quark mixing example 

Yukawa interactions with Higgs fields generate quark masses and mixings. Mixing with 
conventional quarks provides natural decay channels and is expected at some level, since new 
quarks are necessarily unstable |27|] . Suppose first, for simplicity, that mass eigenstates q, x 
arise from the mixing of just one SM quark field q' with a singlet quark field x' of the same 
(unspecified) charge. Then the SM Higgs field H can generate a m'q' L x' R +h..c. mixing term as 
well as the usual mg^g^,+h.c. mass term. A pure singlet mass term Mx' L x' R + h.c. requires an 
isosinglet Higgs field S with vacuum expectation value v$ and coupling (M/vs)Sx' L x' R + h.c; 
this field can also generate a M'x' L q' R term. We then have the 2x2 mass matrix 



( 1 \ 

m m 

v M' M j 



(1) 



where the rows refer to q' L , x' L and the columns refer to q' R ,x' R . This is diagonalized by 
independent rotations of L and R coordinates, giving quark mass eigenstates q and x: 



Ql = q' L cos Ql - x' L sin 9l , qa = q' R cos Or - x' R sin R , (2) 
Xl = q' L sm0L + x' l cos9l , x R = q' R sin R + x' R cos R . (3) 

Since no singlets have yet been discovered, it is natural to assume that the mixing angles 
Ol-i r are small and x is much heavier than q (at least for q = u,d, s, c, b), with m q ~ m, 
m x ~ M 3> m, m', M'. Then g and x are dominated by q' and x' components, respectively, 
with Ol ~ m'/M, R ~ M'/M. Note that SU{2)l gauge couplings relate exclusively to 
and hence are controlled by the left-handed mixing angle 0^ only. 

The heavy mostly-singlet quark x can now decay to q"W and qZ via the couplings 

C xq »w = -A= sin Bl^W^Xl , (4) 
= -ysin ^ cos Ol^l^Z^xl . (5) 

where g is the SU(2) gauge coupling, g z = gj cos Ow and q'[ is the combination of light quarks 
that couple via W to q' L . Since cos 2 0l ~ 1 by assumption, this gives branching fractions in 
the ratio B(x —>■ q"W)/ B(x —>■ qZ) ~ 2 up to phase space factors [^]. Furthermore, if the 
SM Higgs boson is light enough, x can also decay to qH via the coupling 

r 9 m ' - u gsmO L m x _ 

C *« H = ~2M^ qLHXR " 2Mw 1 L R ■ (6) 

Hence the Higgs decay mode too is scaled by sin^, and the three decay branching fractions 
are in the ratios || 

B(x -> q"W) : B(x -> qZ) : B(x -> qH) ~ 2 : 1 : 1 (7) 

up to phase space factors that are close to 1, if m x ^> M w , M z , M H ,m q ,m q /'. These ratios 
can however be altered greatly if this mass ordering does not hold, or if there is large mixing 



In general singlet quarks can mix with all SM quarks of the same charge, requiring a 
more extended formalism. We first consider scenarios with just one new singlet quark. 

(b) One Q — — ^ singlet quark mixing 
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For the case of one charge — | singlet field, mixing with the three SM fields of this 
charge, we denote the mass eigenstate by d, s, b, x where the first three are identified with 
the known quarks (now carrying hitherto unsuspected singlet components) and x is still 
undiscovered. We denote by d' L) s' L , b' L the three orthonormal linear combinations of left 
chiral components that are SU(2) L doublet partners of the known Q = | fields u L ,c L ,t L ; 
the remaining orthonormal combination x' L is an SU(2)^ singlet, and we can write 



( 



d L 



) 



SL 

b L 
\ XL J 



(8) 



V u d V us V u b V ux 

Vcd v cs V cb V cx 
V td V u V tb V tx 

Vod Vos Vob Vox 

Here the 4x4 unitary matrix V generalizes the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
Vckm- The top three rows of V control the SU(2) L gauge couplings of W and Z bosons; 
the first 3 rows and columns of V are precisely Vckm- The submatrix Vckm is generally 
non-unitary. 

The Z couplings to the SU(2) L left-handed doublet and singlet weak eigenstates q' L = 
d' L , s' L , b' L , x' L are given by 

1 



£ = -<7zE <?l fa + - sin 2 9w) rZ,q' L 



(9) 



Hence the FCNC couplings between the mass eigenstates qi = d, s, b, x are 
£fcnc = \g z E ^mYZ,q jL , 



-ij 



<Nu 3 + v*v CJ + v t *v tj 



s a - v oi y oj , 



(10) 

(ii) 



using the unitarity of V . Thus the FCNC coefficients z^- are measures of non-unitarity in 
Vckm- The corresponding FDNC couplings are 

"1 , ,1 



£fdnc = 9z E ^7 m ^m 

i=d,s,b,x 



^za{l - 75) - I sin 2 9 W 



Qi 



(12) 



Thus for the standard d, s, b quarks, mixing with x reduces direct left-handed FDNC by a 
factor (zu — I sin 2 9w)/(^ — § sin 2 9 W ) and leaves right-handed FDNC unchanged. 
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Hence for ra x > M W) M Z) the tree-level widths for CC and FCNC decays to light quarks 



are 



r(x - qi w) 



G F ra\ 
8ttV2 



1 - 



Mir 



rat 



T(x - Qj Z) - Gl '"'''< 



1 - 



2 \ 2 



2M 2 7 \ , „ 



(13) 
(14) 



IGttv^ \ m xJ\ m x / 

If x is heavy enough that all the x — > g^W^ and x ^ qjZ channels are open and all the phase 
space factors are ~ 1, then we can use unitarity to sum over i — u,c,t and j = d,s,b and 
obtain the total CC and FCNC decay widths, 

G F ral 



T(CC) ~ ^Se,|K. 12 



T(FCNC) ~ 



87TV/2 
G F ra\ 



^3\ v oi\ \ v ox\ — 



(1-IK.I 2 ), 



(15) 
(16) 



(17) 



16^ " ° JI 11 " 16^ 
Hence for small mixing (\V ox \ ~ 1) we obtain 

t{cc)/t(fcnc) ~ 2, 

a result proved earlier for two-quark mixing, modulo phase space factors. 

(c) One Q — | singlet quark mixing 

Consider now one Q — | singlet field mixing with the SM fields of the same charge and 
denote the mass eigenstates by u, c, t, x, identifying the first three with the known quarks. 
Let u' L ,c' L ,t' L be the three orthonormal linear combinations of left chiral components that 
form SU(2)l doublets with the known Q = — | fields d^, sl, b F , respectively, while the 
remaining combination x' L is a singlet. We can then write 

^ V u d V us V u b V UO ^ 



U'r c'r t'r x'r — U L C L t L X L 



(18) 



Vcd Vcs Vcb V co 

V td V ts V tb V t0 

y V X d V Xg V X b V XO 

The first three rows of the unitary matrix V control the SU(2)i couplings of W and Z; the 
first three rows and columns of V are precisely Vckm (now generally non-unitary). The Z 



couplings to the SU(2)x left-handed doublet and singlet weak eigenstates q' L = u' L , c' L , t' L , x' L 
are given by 

2 



£ = 9z E <1l {-Ts + 3 sin 2 dw) -fZ^ L 



and the FCNC couplings between mass eigenstates ft = u, c, t, x are 

£fcnc = -Igz^ZzijqiLl^Z^qjL , 

with 

% = v id v; d + v is v* s + v lb v; b = s tJ - v io v* , 



(19) 



(20) 



(21) 



using the unitarity of V. Here again the FCNC coefficients % are direct measures of non 
unitarity in Vckm- The corresponding FDNC couplings are 

r l 2i 

£fdnc = 9z X/ qil^ z v -jzai 1 ~ 75) + g sin 2 ^ ft . (22 

For standard w, c, t quarks, mixing with 2 again reduces left-handed FDNC and leaves right 
handed FDNC unchanged. 



The decay- width formulas are obtained from Eqs.flT3|)- (JT4]), by substituting V xi and 
for Vi x and Zj x . 

(d) One Q = — ~ quark and one Q = | quark mixing 

We here combine the notations of (b) and (c) above, and define 



U L, C £> *L> = C L 5 h, %ul)V 



(23) 



to be three doublet and one singlet Q = | fields, while 



6', 



I' 



d L 



sl 



(24) 
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are the corresponding three doublets (paired with u' L ,c' L ,t' L ) and remaining singlet Q = — | 
fields. Neutral-current couplings of Xd and x u are as in (b) and (c) above. Charged-current 
couplings are defined via the matrix VAV, where A is diagonal with elements 1, 1, 1, 
down the diagonal; VAV generalizes the CKM matrix, its first three rows and columns 



being simply Vckm- For more general mixing parametrizations, see Refs. JTHED 



III. 4x4 MIXING MATRIX 

(a) Experimental constraints 

When extra quarks are mixed in, unitarity constraints no longer apply to the 3x3 
CKM submatrix. Without these constraints, the CKM matrix elements lie in the following 
ranges 



\V\ 



0.9728 - 0.9757 0.218 - 0.224 0.002 - 0.005 
0.180 - 0.228 0.800 - 0.975 0.032 - 0.048 
0.0 - 0.013 0.0 - 0.56 0.0 - 0.9995 



V 



(25) 



/ 



However these numbers were obtained before the evidence for the top quark at Fermilab p9 



The presence of an apparent top quark signal in 6-tagged events at, or even above, the 
predicted SM rate [P9| , pm , strongly suggests a dominant t — > bW decay with \ V t b\ — 1. With 
this extra constraint, all the off-diagonal elements of the 4x4 quark mixing matrix V (or 
V") are necessarily small. One can then generalize the Wolfenstein parameterization 



V us ~\, V ub ~ \ 3 A(p-iri), V cb ~ X 2 A , Vu ~ A 3 A(1 - p - irj) 



by taking for example (for the Q = —\ case V) 



V od ~B(a-ip), V os ~ B(a-ir), V ob ~ B 



(26) 



(27) 



Here the new parameters a, (3, a, r, B are real and B is small (no hierarchy of these 
elements is imposed here). Often it is more convenient to adopt a parameterization in terms 
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of the sines Sj = sin#j of small angles 9i, setting cos#j ~ 1 and neglecting all SiSj terms 
except Sis 2 (see also Ref. [p2|): 



V 



1 

-Si 



-s 3 e tSl + s 1 s 2 —s 2 



1 s 2 s 5 e- iS3 

1 s 4 e~ i52 



\ 



V 



-sse 453 -s 4 e i62 1 



(28) 



In this parameterization, the FCNC coefficients (zij = z*^ are 



Zds = -s 5 s 6 e td3 , z db 



-s A s e e lS2 , z s 



1(82-83) 



(29) 



s 5 e' 



-Ids 



Zbx 



s 4 e 



-M>2 



Similarly, if this parameterization is applied to V in the case of one Q = | singlet quark, we 



have FCNC coefficients (% = 2^) 



183 



-s 5 s 6 e'"*, z ut 



-S4S 6 e l(52 , z ct 



-s 5 e" 



2 fa 



— s 4 e" 



(30) 



-«) 2 



(b) Unitarity constraints 



Unitarity constraints on the 3x3 CKM matrix give linear three-term relations that can 
be expressed graphically as triangle relations in the complex plane; see Fig. |l[ With 4x4 
mixing, they become four-term relations; e.g. for one Q = — | singlet, we have 



VW* + v-v CJ + V*V tj + V*V oj = 5 tv 



(31) 



or again, 



v* d v jd + V*V js + vgv jh + V* X V JX = <%. 



(32) 



For i 7^ j these are expressible as quadrangle conditions in the complex plane. The first three 
terms in each case, however, are precisely the three sides of a triangle if CKM unitarity holds 
(the most discussed example is Eq.([J]]) with % = b,j — d). Thus 4x4 unitarity replaces the 
CKM triangle relations by quadrangle relations. In Eq. (|31~1) the fourth side of the quadrangle 



11 



is V^Voj = —Zij, the FCNC coefficient p|. In Eq.(|3^) the fourth side is V* x V jx , that occurs 
in certain flavor-changing box diagrams (see below). 

In the case of one Q — — ^ singlet quark, the squares of the elements in each row and 
column of the 4x4 unitary matrix V sum to unity. Hence the experimental lower bounds 
on the CKM submatrix elements |28| shown in Eq. (|25|) give constraints: 



1 V ux \ 
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0.08, 


\v cx \ 
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0.57, 


\Vtx\ 
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1.0, 


(33) 
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Also each quadrangle must close, so the exotic fourth side is bounded by the sum of the 
upper limits of the three conventional CKM sides, giving 

r Hi/ I < 0.60, 

(34) 

0.61. 

Finally, when eventually we obtain upper bounds on \V j\ (j = d,s,b) from other data, 
unitarity will imply a lower bound on \V ox \ 2 = 1 — Ej|V^j| 2 , and hence an upper bound on 
£j|^4z| 2 = 1 — \ V ox \ 2 ; this latter bound will apply equally to each \Vi X \ 2 in the summation 
(i = u,c,t). See Section VI (e) below. 

In the case of one Q — | singlet quark, bounds on the CKM submatrix elements of the 
mixing matrix V of Eq. (|i~8|) give analogous constraints: 

\V uo \ < 0.08, \V co \ £ 0.57, \Vto\ < 1.0, 
\V xd \ < 0.15, |t> s | < 0.56, \V xb \ < 1.0, 

( 35 ) 



\v uo \ 


\Vco\ 
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0.44, 




\Vto\ 
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\Vcol 


\Vto\ 


\v xd \ 


\v xs \ 


< 


0.45, 


\v xd \ 


\v xb \ 
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0.03, 


\v xs \ 


\v xb \ 



IV. DIRECT SINGLET-QUARK PRODUCTION CONSTRAINTS 

(a) Z decays 

At e + e~ colliders, xx pairs can be produced directly via their 7 and Z couplings, and 
xq or xq pairs via FCNC. The most stringent bounds at present come from the observed Z 
decay widths, from which it appears that contributions beyond the SM are limited by |H] 
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r z (non-SM) < 15 MeV . (36) 
For the case Q x = — |, the partial widths for decay to one light plus one new quark are 

T(Z -> dx) = V{Z -> dx) = 3r° K \z dx \ 2 F x = (0.66 GeV) \z dx \ 2 F x , (37) 
and for the case Q x = | we have 

r(Z -> uz) = r(Z -> Mi) = 3T° Z K \z ux \ 2 F x = (0.66 GeV) \z ux \ 2 F x , (38) 

where T% = G F M\ / (l2?rv^) = 0.17 GeV , F x = (1 - m 2 x /M 2 z ) 2 (1 + m 2 /2Mf) and K = 
1 + (87r/9)a s (M^) = 1.33 is a QCD factor. For each xq + xq contribution to remain within 
the bound on r^(non-SM) sets m^-dependent constraints on the FCNC coefficients, 



F x \z ix \< 0.11 t = d,s,b, {Q x = -\) (39) 

F x \z jx \ < 0.11 j = u,c, (Q x = §) . (40) 

The partial widths for decays to xx pairs are 

T(Z -> xx) = 24 T° z K [l - 4m 2 /M 2 l 1/2 [<£ (1 + 2m 2 /M§) + (1 - Am 2 x /M 2 z )} (41) 



where 



SV — — \z xx + | sin 2 9w, 9a — — \z xx {Q x — —\), 
9v — \ z xx — | sin 2 0w, 9a = \z xx (Q x — §) • 
In the limit of small singlet-doublet mixing, we have z xx ~ or z xx ~ and hence gy ~ 
— Q x sin 2 9 W , 9a — 0. In this limit the upper bound T(Z — > xx) < r^(non-SM) gives 

m x > 42 GeV (Q x = -§) , (42) 
> 45 GeV (Q x = |) . (43) 



m 



X ~ 



Figure |2| shows the corresponding xx contribution to r^(non-SM) versus m x for Q x = — | 
and |. Direct searches at LEP for typical heavy quark signals (t — > 6' — > cW - ). 

based simply on event shapes, set early limits m t > 44.5 GeV and my > 45.2 GeV 



corresponding to upper limits YiZ — > b'b',tt) < 20 — 30 MeV. Applying these limits to 
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singlet quarks gives a weaker result than Eq.(f42|) and about the same as Eq. (f43|) ; some 
improvements could presumably be achieved with present much higher luminosities. 

(b) Hadroproduction 

At hadron colliders, xx pairs can be produced via QCD interactions exactly like SM 
quark pairs. Their x — *> q'W CC decays into lighter quarks give signals rather similar to 
the t — ► bW signals that have been looked for in top-quark searches p9j[30l , [33f , although the 
details may differ; they also have new decays into qZ and/or qH. We briefly discuss some 
examples. 

(i) For a heavy Q — | singlet Xt that mixes mostly with t and has Mw < Tn Xt < Mw + m t , 
the dominant decay mode is Xt — > bW while xt — > tZ, tH are kinematically forbidden. Hence 
the x t x t signals look exactly like it signals, including the presence of taggable 6-jets in the 
final state. Lower bounds on mt such as the DO result |K| mt > 131 GeV apply also to m Xt . 
Recently published evidence for it production |29] could in principle be interpreted as XtXt 



production, but electroweak radiative corrections |34| already indicate a top mass near the 
observed value, making it production the most likely interpretation. However, if there is an 
excess of top-type events above the SM rate this could be due to x t x t production in 



addition to it production 

(ii) A Q = — | singlet Xb that mixes mostly with b and has M z < m Xb < m t + M w would 
decay dominantly via x& — > bZ, bH with the tW mode suppressed, escaping the usual top 
searches but offering new Z and H signals. If the latter are suppressed (e.g., if m# > m Xb ), 



early CDF limits on the remaining Z signals imply a bound m Xb > 85 GeV [|T^]. This 



scenario gains fresh interest |L6| from hints of possible excess tagged Z plus four jet events 
at the Tevatron p9| . 



To be quantitative about signal expectations with 6-tagging, let us consider x fe and t to 
be degenerate (m Xb ~ m t ) for simplicity, so that they are produced equally. If Xb is lighter 
than this, the singlet signal rates will be correspondingly higher. For the singlet decay we 
consider two extreme scenarios: (A) tuh > rn Xb > Mz, so that T(xb — > bZ) ^> Y{xb — > bH) 
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and (B) ran — Mz with small b-Xb mixing, so that T(xb — > bZ) ~ T(xb — ■» 6-ff"). For 6-tagging 
efficiency, we assume e& = 0.2 to be the probability for tagging a single fe-quark; then the 
probability for tagging a bb event is 1 — (1 — €b) 2 = 0.36 and the probability for tagging a 
bbbb event is 1 — (1 — e^) 4 = 0.59. Our discussion is simplified in that we neglect fake tags, 
assume 6-tags are uncorrelated, and assume 100% acceptance. Then the probabilites for 
different final state configurations including 6-tagging are 

channel probability with tag 

bbWW -> bb{£u){jj) 0.29 x 0.36 = 0.104 

bbZZ -> bb{££){jj, bb) 0.094 x 0.41 = 0.039 

bbZH -> bb(££)(bb) 0.067 x 0.59 = 0.040 

summing over £ = e, fi channels. The first numerical factor on the right is the branching 
fraction and the second factor is the 6-tagging probability. Thus the leptonic W/Z event 
ratios in our two m Xb ~ m t scenarios (A) and (B) are 

N(tt -> W iv + 4j with tag)/N(xx -> Z« + 4j with tag) ~ 2.7(A) or 3.5(5). (44) 

In contrast, the QCD electroweak background ratio is B5f 



N(QCD -> W fe + 4j with tag)/A^(gCD -> Z« + 4j with tag) ~ 10 - 14. (45) 

(iii) A Q = — | singlet quark x d mixing mostly with d would decay by x — > uW, dZ, dH in 
the ratios 2:1:1 modulo phase space factors. Thus for m x 3> Mw, Mz, Mh the top-like 
signals would be reduced roughly by a factor 2 for single-lepton channels and by a factor 4 
for dilepton signals, compared to a top quark of the same mass; however, for smaller m x the 
reduction is generally less, and in the window Mw < rn x < Mz, Mh there is no reduction. 
But there is now no 6-quark to tag. Examination of earlier top-quark searches without a b- 
tag P0y53|) scaling down the top-quark expectations by some factor between 1 and 4, shows 



that the range M w < m x < M z , M H is definitely excluded, and probably some adjacent 
ranges of m x too, but more cannot be said without detailed analysis. Similar conclusions 
apply to x s singlets mixing mostly with s and to charge | singlets x u or x c mixing mostly 
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with u or c; there are small differences between these cases, such as the lepton spectrum flO 
and the taggability of c-quarks, but they do not change the overall conclusion, 
(iv) Decays outside the detector. The distinctive x-quark signals will be lost if x decays 
outside the detector. If we assume typical Lorentz factors j3 , y ~ 2 and require that the mean 
decay distance lp = P^c/T due to any single x — > qtW, qjZ decay mode exceeds one metre, 
Eqs.dD-dH give 



|VU < 1.2 x 10" 



for Q = — \ , and similarly 



ivy z 1.2 x lo- 



200 GeV 



\V oj \ < 0.9 x 10" 



200 GeV 



m. 



(46) 



200 GeV 



V oj \ < 0.9 x 10" 



200 GeV 



(47) 



for Q — |. If these conditions hold for all light quark flavors i,j , then singlet decay signals 
at hadron colliders will be greatly suppressed. [The conditions are somewhat weaker for 
m x < M W ,M Z }. 

(c) Leptoproduction 

At ep colliders, singlet quarks can be produced by the same jg fusion processes as SM 
quarks; Zg fusion is also possible (at reduced rates due to reduced Zxx couplings) but Wg 
fusion is only possible via mixing. However, the reach of the HERA collider for new quark 



detection is much less that that of the Tevatron []36| so this is not a promising avenue for 
singlet discovery. 

(d) Summary 

The LEP mass bounds Eqs.(|42|)- (|43D are virtually unconditional. Hadroproduction 
bounds are much stronger in particular cases [e.g. m x > 85 GeV (131 GeV) if the decays 
x — > qZ (x — > qW) dominate completely] , but assume implicitly that the mixing with lighter 
quarks is not so extremely weak that x decays outside the detector [typified by off-diagonal 



4th row and column mixing-matrix elements all being < 10 8 [(200 GeV)/ 



m x )2 
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V. NEUTRAL MESON- ANTIMESON OSCILLATIONS 



The existence of new singlet quarks can affect neutral meson-antimeson oscillations in two 
different ways, through FCNC tree-level Z exchange and through box diagrams, illustrated in 
Fig. [5] for the B®-B® case. In the tree diagram of Fig. |3|(a), the effects are due to an additional 
Q = — | singlet that generates the FCNC couplings of Eqs. (|ToD -(|Tll); more generally, such 
FCNC effects of Q = — | singlets occur also for K°-K° and B Q S -B® oscillations, whereas 
analogous effects from Q = | singlets give D°-D° oscillations. 

In the box diagram of Fig. |3|(b), the effects come from an additional Q = | heavy quark 
option in the loop, along with corresponding reductions in the original three generation 
couplings. Similar effects are present in K°-K° and B° s -B° s oscillations, while analogous 
effects from Q = — | singlets occcur in the D°-D° case. The association of singlets with 
their mixing effects is summarized in Table |. 

(a) Z-exchange contributions 

We first analyze the Z-exchange FCNC effects, which are potentially the most interesting. 

For B®-B® oscillations, the contribution is 

V2G F m B f^B B r] B 2 
l<H = 3 1**1 , (48) 

where fs is the B^ decay constant, Bb is the bag factor (B — 1 is the vacuum satu- 
ration approximation) and t]b ~ 0.55 is a QCD factor. (We assume the QCD correc- 
tion is the same for both the Z-exhange contributions and for the box diagram contri- 
butions described below.) The analogous expressions for K°-K°, D°-D°, B° s -B® oscilla- 
tions involve Re(z^ s ), z^ c , z 2 sb) respectively. Actually this FCNC process contributes co- 
herently with the SM box diagrams. We shall here assume very conservatively that the 
singlet-quark Z-exchange contributions do not exceed the measured values. Then from 
the measurements \5m\ K = (3.51 ± 0.02) x lfr 12 MeV, \5m\ D < 1.3 x l(T 10 MeV, 
\Sm\ Bd = (3.4 ± 0.4) x 1(T 10 MeV we obtain 

\Re(z 2 ds )\ = \{Rez ds f - (Imz ds ) 2 \ < 9 x 10" 8 , (49) 
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(50) 



\z db \ = \V od \ \V ob \ < 8 x 10- 4 , (51) 
where the < symbol reflects some uncertainties in the factors /, B, rj. Similar bounds on 



z ds and z db are given in Ref. |Jl8| , p2| , |25| l) and on z uc in Ref. p6[ . For |<5m| Ba there is only an 
experimental upper limit [^] and hence no bound on \z s i,\ = \V os \ \V b\. We take the D and 
B decay constants from Narison f38|: fn = 1.37 f n , = 1-49/,,-, with /„- = 0.131 GeV and 
fx = 0.160 GeV, and set B — 1 and rj = 0.55 in all cases. Taking the lower bound B = | 
instead would raise the limits above by a factor \f?>. 

(b) New box diagram contributions 

In the case of box diagram contributions, the constraints on the mixing are rather differ- 
ent. First consider the case of B^-B® box diagrams with an additional Q x = | contribution, 
with m x ~ m t approximate degeneracy. Then the SM formula is 



G\ B f 2 B m B ri B 
6tt 2 



\Sm\ S M = — J "7 a ' m I^^Ickm l J sl > (52) 



where 1b is a box-integral factor (see e.g. Ref. |39|), and the effect of adding an extra singlet 
is to replace the CKM factor by 



v xd v: b + v td v t ; 



2 2 



V ud V* b + V cd V* b \ CKM (53) 



using unitarity. However, \V td V t * b \cKM = \V ud V* b + V^V^Ickm, so the prediction for \8m\ is 
effectively unchanged in this x, t mass-degenerate limit. Only if x is much heavier than t 
can significant changes arise. Similar conclusions apply to K°-K° oscillations. 

The Q — — g box diagram contributions to D°-D° mixing are potentially more interest- 
ing, because d, s and b are relatively light compared to the allowed mass scale for x. Here 
the x contributions may be dominant (depending on the size of the mixing) and given by 



\Sm\ D = — \VcxV ux \ \I D \ , (54) 

07T Z 



where — m 2 , for m x ~ 200 GeV, giving 



\V CX \ |K,|<0.7xlO- 2 f^-^) . (55) 



A similar bound is noted in Ref. [|D| for mixing a fourth-generation V quark, that is essen- 
tially equivalent to singlet mixing in this context. It is expected that a future sample of 
10 8 reconstructed .D's would have a factor 20 improvement in sensitivity to Smp, and would 
consequently give a factor ~ 4 — 5 more sensitivity to the above mixing. Note that SM 
short- and long-distance contributions are far below this sensitivity f41~l , |42| . 



The parameter ex , that describes CP-violation in K° — K° oscillations, also receives 
tree-level Z-exchange contributions from Q = — | singlet mixing: 

i i G F m K B K f\ 2 

\e K \ = ul6mKl \Im(z ds ) |. (56) 
Requiring |e^| < |e^| ex . p = 2.27 x 1CT 3 gives the bound [p^ .^5[ 



\Im(z ds ) 2 \ < 6 x 10- 10 , \Re{z ds ) Im{z ds )\ < 3 x 1(T 10 . (57) 
Combined with Eq. (|49|) , this gives 

\z ds \ = \V od \ \V os \ < 3 x 10" 4 . (58) 



VI. FCNC DECAYS AND FDNC EFFECTS 



(a) Kl — > decay 

Experimental measurements on FCNC decays imply constraints on the FCNC Z cou- 



plings and hence on singlet-quark mixing parameters [p|,^|, p!8| , p2| , p3| . For example, Kl — > 
fi + fi~ has a Z-mediated diagram if a Q — — | singlet x mixes with d and s, contributing the 
decay width 



r(K L 



V V )z 



2G I F ff ( m K m 



8tt 



^[l~4ml/m 2 K ] 1 / 2 [(--sm 2 e w ) 2 + (sm 2 e w ) 2 } \z ds \ 2 . (59) 



After subtracting the contribution for the 77 intermediate state (an imaginary decay am- 



plitude), the latest Brookhaven results |13] indicate an upper limit B rea i < 5.6 x 10 (90% 
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CL) on the contribution to the branching fraction from the real part of the decay amplitude. 
Using this to bound the contribution from Re(z ds ) we obtain 

\Re(z ds )\ <0.64x 10~ 5 . (60) 

The combined bound on \z ds \ remains unchanged. 

(b) B°, D° -> fi+fx- decays 

Analogous formulas describe the tree-level contributions to D° — > fi + p~ and B° — > 
decays (without the factor 2 because D° and B° are not pure CP = — 1 states). Requiring 
that the Z-exchange contributions are within the experimental limits B(D° — > p + < 
1.1 x 10" 5 and B(B° -> < 5.9 x 10~ 6 B7| , gives the constraints 



\z 



= \V uo \ \V m \ < 0.20 , (61) 
za\ = \Vod\ \V ob \ < 0.04 , (62) 



much weaker than the oscillation bounds Eqs.(|5l3|)-(|5"1|). 

(c) B, D -> Xg+T; decays 

The rare decays B — > X£ + £ _ occur at tree level, via FCNC couplings and z sb which 
give 

r( J B^£+z/A) 2 |Kb|+p|K6| 

where p ~ 0.5 is a phase space factor; p = 1 — 8r 2 + 8r 6 — r 8 — 24r 4 ln(r) with r = m c /m b = 
0.316 ± 0.013. Hence the experimental limit B(B — > Xp + p~) < 5.0 x 10~ 5 gives the 
constraints p2"U?3 | 



\z db \ = \V ob \ \V od \ < 0.04 x \V cb \ < 2 x 10~ 3 , (64) 

\z sb \ = \V ob \ \V os \ < 0.04 x \V cb \ < 2 x 10- 3 . (65) 

The first bound is competitive with that from B d — B d oscillations in Eq.flSTp. Upper limits 
have recently been given for some D — > p + p~ -l-hadrons branching fractions [Q, suggesting 
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an inclusive upper limit of order (1 — 2) x 1CT 3 (although no explicit value is quoted); such 
a limit would however only give \z uc \ < 0.2 — 0.3, possibly competitive with Eq. ( |6"ID but 
much weaker than Eq. (|50|) . 

(d) B — > s(d)j decays 



The rare decays B — > s(d)7 have also been considered f24l , f45"l| . In the SM they go via 
VT-loop diagrams; adding a down-type singlet quark introduces new Z-loop diagrams, using 
the FCNC couplings (if-loops are usually negligible). These can be incorporated into 
the conventional analyses by adding their contributions into the coefficients of the effective 
operators of the magnetic and chromomagnetic moment couplings fjp and as described 
in Appendix A. The ratio of T(b — > qj) (where q = d, s) to the inclusive semileptonic decay 
width is then given by 



r(6- 


-> <n) 


6a 


\v t ;v tb \ 2 


r(6 - 


■> cev) 


iipX 


\Vcb\ 2 



\c 7 (m b )\' 



(66) 



where a is the electromagnetic coupling and 
c 7 (m b ) 



~a s {M w )~ 


16/23 | 


_ a s (m b ) _ 





c 7 (M w ) - -c 8 (M w ) 



'a s {M ] 



Wj 



' a s {m b ) 

M M w) 



2/23' 



(67) 



i= i \ a> s (m b ) J 

The Wilson coefficients c 7 and eg, the coefficients hi, and the exponents a* from the 8 x 
8 anomalous dimension matrix ^6[ are given in Appendix A. The phase-space factor p 
is defined below Eq.(^3]) and the QCD correction factor A for the semileptonic process 
1 - 2/(r,0,0)a a (m 6 )/7r with f(r, 0, 0) = 2.41 0. We remark that the FCNC 



is A 



diagrams include not only Z-loops but also tree-level Z-exchanges between the b-quark and 
the spectator antiquark in a decaying B-meson, not commented upon in previous literature. 
However, these Z-exchanges are suppressed relative to Z- loops by factors fs/n^B ~ 1/25 in 
decay amplitudes [jUjfl , so we do not pursue them here. 

In the SM one expects the ratio B(b — ► d^)/B(b — > 57) w |Vtd/Vf S | 2 , since the QCD 
corrections largely cancel out. The additional FCNC terms are proportional to z qb /(V tb V t *) 
in each case ( 



d, s), and it has been shown that |24 
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These limits permit singlet quarks to have greater impact on the b — > d'-f rate (e.g. if 
Zdb ~ z sb)- On the other hand, one expects from the general decoupling theorem [|9[ that 
Zdb is much smaller than z sb . 

An up-type singlet quark can also be considered. Its contribution is the same as from a 
standard fourth generation, giving 



|2 



(69) 



T(b — > cev) irpX \ \V c b\ 

where the matching conditions for the relevant Wilson coefficients are again given in Ap- 
pendix A. The major contributions to B{b — > 57) are now the t- and x-quark loop terms 
in Eq . (|69D . Notice that T(b — > 57) is the same as in the SM when m x = m t , for the same 
reason as in B®-B® and K°-K° oscillations above. But if m x deviates significantly from 
m t , an enhancement or suppression relative to the SM can be expected (as with a fourth 



generation |50|). 

Figure £| shows the b — ► S7 rate versus m x with various values of |V^*t^. 6 |, for the SM 
plus one up-type singlet quark. We have assumed here that the phase of V* s V xb is the same 
as that of V* s V tb within a sign, so that deviations from the SM are maximized. We note 
incidentally that the unitarity constraint on IV^I |Krb| helps to guarantee that B(b — > dj) 
with a w-type singlet remains close to the SM. 

(e) Z decays 

We turn now to FDNC effects. At tree level, introducing mixing with a singlet quark x 
simply reduces the left-handed coupling of a conventional quark i by a factor 1 — | V^| 2 /(l — 
! sin 2 6^1/ ) (for charge Q x = Qi = — |), or by a factor 1 — |Vi G | 2 /(l — ^sin 2 8 w ) (for charge 
Qx — Qi — §), leaving right-handed couplings unchanged; see Eqs. ([12]) , (|20|) . We shall 
neglect singlet-mixing effects at one-loop level, where they are small corrections to small 
corrections. 

The Z partial decay widths, branching fractions and asymmetry measurements directly 
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probe the FDNC Zqq couplings. Z — > bb decay is an interesting case to consider, since there 



is at present some discrepancy between the LEP data ]3T] and the SM prediction for the 
ratio Rb = T(Z — ► bb)/T(Z — > hadrons): 

i^(LEP) = 0.2202 ± 0.0020 , R b (SM) = 0.2156 ± 0.0004 . (70) 

Since b — x mixing reduces the Z — > bb coupling, it would make the discrepancy worse. The 
decay width has the form 

nz - ib) = y^iMlp j>(^ )2 + , (7i) 

where (3 is the CM velocity and g\ and g\ are the axial and vector Zbb couplings, so down- 
type singlet mixing dilutes the tree-level contribution by a factor ps (1 — 2.4|V ;,| 2 ). It is 
inadvisable to derive a limit on V Q b from this result alone, however, since the SM itself is on 
the verge of being excluded. Many models with down-type singlets also give corrections to 



Z — ► bb from mixing Z with a new Z'; these too are typically negative pT |. 

A global comparison of all FDNC effects with the latest LEP and SLC data leads to the 
following constraints (see final paper of Ref. [pl"H): 

\V od \ 2 < 0.0023 , \V os \ 2 < 0.0036 , \V ob \ 2 < 0.0020 , 
|K | 2 < 0.0024 , \V co \ 2 < 0.0042 , 

assuming at most one singlet quark mixes with each conventional quark. From these num- 
bers, unitarity of V then gives 

\V ox \ > 0.996, \V qx \ < 0.089, {q = u,c,t). (72) 

(f) Other FDNC effects 

Singlet mixing could also change FDNC effects in neutrino scattering and atomic parity- 



violation measurements [20], but there appear to be no useful constraints from this quarter. 



23 



VII. CP ASYMMETRIES 



The amount of CP violation in the SM is measured by the size of the unitarity triangle 
in Fig. |l[ How this CP violation shows up in decays is determined by the angles of the 
unitarity triangle(s), which appear as CP asymmetries in decays to CP eigenstates. The 
angles 

a ( v °* v *\ - ( Y^k\ ™ 

"^nwaJ- a=arg o^o)' (73) 

that characterize CP violation, are directly measurable in Bd decays with b — > c and b — > u 
respectively. The prototype processes for measuring (3 and a are Bd — > tpKs and B& — > 



+rrr respectively. [The angle 7 = &rg(—V u dV* b /V c dV* b ) can be measured in the decay 



7T 1 7T 



B s — > pi^5, which will prove much harder at a B factory because of the small branching 
fraction and the possible contamination from penguin contributions.] Present information 
on the third generation couplings does not tell us much about the asymmetries. Future 
improved measurements of the CKM mixing angles will pin down the SM prediction more 
precisely. We find the biggest uncertainty in the SM asymmetries stems from the uncertainty 
in V u b, a quantity ripe for better measurement at a 5-factory. 

We assume as usual that the asymmetries are dominated by the interference between 
two amplitudes, one of which is given by B d -B d oscillations with Ti2 << Myi- The time- 
dependent CP asymmetry in the decay of a B® or Tf d into some final CP eigenstate / is 

vm) - /) - r(gw -* /) = _ Im ^ } t) 

r(B2(t)-/) + r(lto-/) 

where 5m is the (positive) difference in meson masses, the mesons states evolve from flavor 
eigenstates B® and B° d at a time t = 0, and Im X(Bd — > /) is the time-independent asym- 
metry. The quantity Im \{B d — > /) is — sin 2(3 and sin 2a for / = ipKs and / = 7r + 7r~ 
respectively in the SM (we neglect possible penguin diagrams in the decay B d — ► 7r + 7r _ . 

We consider the allowed range for the Wolfenstein parameterization involving p and 77 
recently given in Ref. |5T] . The angles a and /3 are easily related to p and 77 through the 
unitarity triangle in Fig. |l| 
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sin 2a 
sin 2/3 



2rj ( V 2 + p(p-l)) 
(r/ 2 + (1 - p) 2 ) (r/ 2 + p 2 ) ' 
277(1 -p) 



(75) 
(76) 



p 2 + (l-p) 2 ' 

In the presence of d-type singlet quarks the unitarity triangle becomes a quadrangle 
as described in Section III, and the CP asymmetries in B decays are altered from SM 
expectations. The deviations occur in two ways. 

(1) The angles j3 and a no longer have SM values, because the revised unitarity constraint 
yields different allowed ranges and more general phases for the CKM elements. 

(2) There is an additional B d — B d oscillation contribution from tree- level Z- mediated graphs. 
The asymmetry expressions are modified to 



Im X(B d -> ^K s ) = - sin (2/5 + argA M ) 
Im X(B d — > 7r + 7r~) = sin (2a + argA M ) , 



(77) 
(78) 



where 24 



A bd 



1 + r d e M 
4vrM^ sin 2 6 



9 bd = arg 



aI B {x t ) 



z bd 



VuV t * b 
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z bd 



VuV t * b 



and Isi^t = m t/^w) * s the box integral(see e.g. Ref. |39| 



lB{x t ) 



1 



w 



art 1 + 
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6xf 



Inxt 



(79) 
(80) 

(81) 



(82) 



x t (l-x t ) 2 ) (1-XtY 

The contribution of z d b to the unitarity quadrangle can be described by a magnitude 
and a phase 9m (relative to V td y* b ). This phase can take any value between and 27r, but 
the magnitude must be consistent with closure of the quadrangle. In Fig. [| we show the 
asymmetry for the decay B d — > ipK$ for different values of the parameters [E4I 



Zbd 

VuV t t 



9 



bd 



(83) 
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with the CKM angles, the top mass and the B lifetime fixed at their central values: \V u d\ = 
0.9743, \V cd \ = 0.204, \V ub \ = 0.0035, |V c6 | = 0.40, m t = 174 GeV and r B = 1.50 ps. We 
take \J Bfg = 195 MeV as we did previously, and use the next-to-leading order value for the 
QCD correction tjb = 0.55 ||52|| . By taking the coherent sum of the contributions to B^-B® 
mixing from Eqs. ( ^8|) and (|5^), and the the mixing parameter Xd = Sm/T = 0.71, one can 
determine |V^V^|. The shaded band in Fig. ||] indicates the allowed range in the SM for the 
asymmetry Im X(Bd — > ipKs)- 

The quantity 5d can be quite large, as indicated by Eq. (|68|), but Fig. [| shows big effects 
even with much smaller Sd- One notices that the CP asymmetry Im X(Bd —>■ ipKs) is 
negative in the SM, but with sufficiently large 5d one can obtain positive values p2LE4|. The 



effect of singlet quarks on CP asymmetries can therefore be dramatic [22,24 



The CP asymmetry Im X(Bd — > 7r + 7r~) is shown in Fig. |^ for various values of 5d and 6m- 
Here the SM expectation covers the entire range, so merely measuring the sign of the CP 
asymmetry could not upset the SM. But given well-determined CKM elements, deviations 
from SM predictions could be significant and could provide evidence for singlet quarks. 



VIII. GUT SOURCES OF SINGLET QUARKS 



(a) Generalities 

GUT models provide arguments for the existence of particles with exotic quantum num- 
bers, but also impose restrictions upon them. In this Section, we explore the constraints on 
singlet-quark models implied by coupling-constant unification and perturbativity. Most of 
the examples we consider are supersymmetric models, and one must bear in mind that these 
models have extra contributions to the processes described above, so that the constraints 
obtained can be affected. 

Singlet quarks considered alone do not introduce gauge (or gravitational) anomalies, but 
they spoil the successful gauge coupling unification of the minimal supersymmetric model 
(MSSM) if the singlets are below the GUT scale, since they change the running of the 
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SU(3) and U(l) couplings but not the SU{2) coupling. For down- type singlets, this can be 
remedied by adding more fermions to fill out the 5 and 5* representations of SU(5) or the 
10 of 50(10); see the examples below. 

For up-type singlets, however, it is harder to find a consistent scenario, if one believes 
that gauge coupling unification is due to a GUT symmetry and therefore wants to retain the 
desert between the GUT scale and the scale of the exotic fermions. The model of Barbieri 
and Hall f| postulates that singlet quarks arise as supersymmetric partners (gauginos) of 
gauge bosons from a unification group that assigns a fourth color to leptons, so here the 
singlet with the right quantum numbers to mix with the top quark is not a matter fermion at 
all. We can introduce top-like singlets as matter fermions by assigning them to the adjoint 
representation of the GUT group. The smallest suitable representation of matter fermions 
is then the 45 of SO (10), or the 78 of E$. But these representations are too large; they 
destroy the asymptotic freedom of the strong coupling, and contain extra doublet quarks 
besides. Alternatively, in the context of SU(5), we can introduce one up-type singlet quark 
by adding one extra light 10 and one 10* representation; these bring one extra vector-singlet 
lepton plus a vector-doublet of quarks too, and restore gauge unification with 6 3 = at one- 
loop level. Two-loop effects become large, however, and large threshold corrections must 
be invoked to restore gauge coupling unification. Apart from this 10 + 10* scenario, there 
appears to be no simple way to arrive at a low-energy model with up-type singlet quarks 
from a desert GUT model. 

(b) One-loop evolution equations 

The evolution equations for the gauge couplings at one loop can be written 



dt 16tt 2 ' dt L 1 1 2tt 
with t = ln(n/Mz) the logarithmic scale and = gf/(47r). The SM particle content alone 
gives b\ = 4^j, b 2 = — 3|, 63 = —7. It is well known that this does not lead to gauge coupling 
unification; given a^ 1 , af 1 evolves too fast compared to a^ 1 . Simply adding singlet quarks 
makes things worse, however; a^ 1 is unchanged, aj -1 evolves faster and a^ 1 evolves more 
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slowly. Some different particle content is needed. 

The MSSM, with additional supersymmetric particle content and two Higgs doublets, 
does give successful gauge coupling unification with the beta functions 

h = 2n G + ~n H , (85) 
5 

b 2 = 2n G + n H - 6 , (86) 

6 3 = 2n G - 9 , (87) 

where n G = 3 is the number of light generations of matter and nn = 1 is the number of pairs 

of light Higgs doublets. In the presence of n Xu up-type and n Xd down-type singlet quarks, 
the beta functions are modified to 

3 2 8 

61 = 2n G + -n H + -n Xd + -n Xu , (88) 

b 2 = 2n G + n H -6, (89) 

63 = 2n G - 9 + n Xi + n Xu . (90) 

These singlet quark contributions upset the MSSM unification as shown in Fig. [^. However, 
unification can be restored by adding exotic fermions to fill out one or more representations 
of SU(5). For example, in an E e model the basic 27 representation has the decompositions 
16 + 10 + 1 in the 50(10) subgroup, which in turn are (10 + 5* + 1) + (5 + 5*') + 1 in 
SU(5). If the full 27 of fermions is light, the down-type quarks are supplemented by colorless 
doublets and singlets, giving 

3 

h = 3n G + -n H , (91) 
b 2 = 3n G + n H -6, (92) 
63 = 3n G - 9 . (93) 

where n G is now the number of light generations of E$ matter, assuming that the light Higgses 
are external to the 27 representations. Thus all bi are shifted by the same amount n G from 
the MSSM case; all values of a~ l (t) are shifted down by the same amount —bitn G /(2n) and 
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unification is preserved. Note incidentally that with the usual three generations we now 
have 63 = and the strong coupling ceases to run. 

However, one usually assumes instead that the pair of Higgs doublets comes from the 
27, in which case the beta functions are 

61 = 3n G , (94) 
b 2 = 3n G - 6 , (95) 
b 3 = 3n G - 9 , (96) 

and the gauge coupling unification is again problematic. One solution would be to get back 
to the previous successful beta functions by adding two new particles with the quantum 
numbers of two Higgs doublets. Alternatively, we might have hoped that the two-loop 
contributions could rescue gauge coupling unification, since the gauge couplings are all larger 
than in the MSSM, making two-loop contributions more important. Unfortunately the sign 
of the two-loop term pushes the SU (3) coupling further away from the electroweak couplings. 
Nevertheless, this example shows that two-loop contributions could be important. 

(c) Two-loop evolution equations 



At two-loop level, the evolution equations become 

dgi 9i 



dt 16tt 2 
where the one-loop beta functions are 



1 3 

U 2 + ^DA 2 



(97) 



3 1 3 

h = -n w + -(n 5 * + n 5 */ + n 5 ) + -n H , (98) 
II 

3 1 

h = ~n w + -(n 5 . + n 5 */ + n 5 ) + n H - 6 , (99) 
3 1 

63 = 2^10 + -(n 5 * + n 5 */ + n 5 ) - 9 , (100) 

for an arbitrary number of copies (n w , n^*, 715*', 715) in the 10, 5*, 5*', and 5 representa- 
tions of SU(5), and n H light pairs of Higgs doublets (from a split representation) . The 
two-loop coefficients bij are listed in Appendix B. The model-dependent contributions from 
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the Yukawa couplings at two-loop order have been neglected. In the absence of split repre- 
sentations the entries are related by simple SU(3) and SU(2) group factors. The second 5*' 
representation contains the multiplet (H, x^) which may come from either the 10 or the 16 
representation of 5*0(10). 

For rank 5 Eq models there is an extra U(l) that enters into the gauge coupling evolution 
equations at the two-loop level. There is a one-parameter family of extra t/(l)'s orthogonal 
to U(1)y- Three popular models [ 53| are characterized as follows: 



(1) the SO (10) singlet fermion is inert, with respect to the extra £7(1). 

(2) the SU(5) singlet fermion in the 16 of 5*0(10) is inert, and 

(3) the 5* and 5*' have exactly the same quantum numbers. 

We label the U(l) quantum number of these models by Y', Y", and Y'" respectively, and list 

the quantum numbers for the full 27 of E§ in Table 0. Notice that for the first model the 

16 of 50(10) decomposes as 10 + 5*' + 1. The first two models could actually arise from 

an 50(10) theory, since Tr Y vanishes across each 50(10) multiplet, while the third model 

is distinctively E 6 (as are all the rest of the rank 5 models). In the last model it is natural 

for the entire 27 to be light. [We note that the extra abelian groups are often referred to 

by the notation U(l) v , U(l) x and U(l)^ [Q; the models (2) and (3) considered here then 

correspond to the extra group being U(l) x and U(l) v respectively, while the model (1) is a 

linear combination of U(l) x and U(l)^.] For these three models one obtains, in addition to 

the RGE coefficients bi and b^ already listed, the one-loop coefficients 

,/ 1 1/ 9 5 13 /ini . 

b p = -nio + -{n 5 * + n 5 ) + -n B .# + -n N + —n H , (101) 

19 1 5 2 

K = -n w + -n 5 * + -(n 5 .# + n 5 ) + -n v c + -n H , (102) 
y 4 8 2 8 ^5 

... 2 1 . .45. ,17 . . 

b P = 3^10 + Y2^ 5 * + nb *'> + 3 ns + l2^ nN + Uv V + 30 nH ' ^ 103 ^ 

and the two-loop contributions are listed in the Appendix B. We use the subscript p (for 
"prime") to distinguish the U(l)', U(l)" and U(l)"' gauge couplings and their RGE coeffi- 
cients. Here and n u <^ are the number of light singlets with the quantum numbers given 
in Table ||. In a general model with more than one U(l) factor, one must account for the 
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mixing between the ?7(1) 's in the renormalization group equations [j54 ]. This complication 
does not arise if one considers only unification trajectories. In any event the practical effect 
of the extra U(l) on the evolution is small. The above equations have been derived for an 
arbitrary number of different representations of SU(5), but split representations in the 10 
of S'O(IO) have been allowed for (10 — > 5 + 5* — > H + H), since they may be needed to 
achieve gauge coupling unification. The beta functions for U(l)' and U(l)" are related by 
ris* <-> ri5*' and un <-» w^. One can also consider the continuous family of rank 5 E§ models 
that include the three above, but as far as gauge coupling unifications is concerned they 
offer no new features. 

An Eq model with three light generations would have nio = n^* = n u ^ = n 5 = n 5 */ = 
tin = 3 (from the usual decomposition of the 27 representation). If only complete 27 
multiplets of E 6 occur, the above coefficients become universal, namely 



and 



b 1 . = b"- = b'" 

13 l 3 *J 



h 


= K = 


K = 


h 


= 3n G 


-6, 


bs 


= 3n G 


-9, 


U 


1 3 


8 


1 


3 3 


8 


1 


1 21 


8 


K 1 


1 3 
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3n 



G 



(104) 
(105) 
(106) 



n G 



-24 



-54 



(107) 



(d) Specific scenarios 

These results allow us to examine gauge unification with different numbers of light gen- 
erations of exotic matter, with or without a pair of light Higgses from a split representation. 
We recall that the MSSM has uiq = n^* = 3, = n^*i = with nn = 1 (the distinction 
between the 5* and 5*' is immaterial as far as SU(5) is concerned), and the Higgs contribu- 
tion is vital for successful unification. 
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(i) Extensions of the MSSM. Adding just one or two generations of SO (10) 10-plet matter 
(72.5 = n 5 *i = 1 and keeping i%h = 1) yields successful gauge coupling unification as shown 
in Fig. |8|. With three generations, however, the two-loop corrections threaten to spoil uni- 
fication; they also make increase toward the GUT scale (although 63 = at one-loop). 
On the other hand, one expects the low-energy threshold corrections to be more significant 
in this case [^], and ultimately the success of unification depends on the details of the 



low-energy spectrum. [Ref. [|55[] assumes that the 577(5) multiplets are degenerate at the 
GUT scale.] 

(ii) Eq based models. Here one needs some split representation since otherwise the elec- 
troweak couplings do not run fast enough for successful unification. (Some attempts at E 6 
phenomenology have assumed that the Higgs pair comes from a complete light representa- 
tion.) The two-loop contributions do not help since they tend to slow the running of the 
strong coupling constant, or even make it grow in the case where the one-loop beta function 
63 is exactly zero (as happens for three generations of light E$ matter). Although asymptotic 
freedom is lost above the exotic fermion mass scale, this is not necessarily a problem for 
gauge coupling unification as long as the two-loop effects do not make a 3 nonperturbative 
below the GUT scale. One has gauge unification at the same scale as in the MSSM (neglect- 
ing threshold correction), with unification coupling a3(M exot i C ) still perturbative, though 
significantly larger than in the MSSM. 

(iii) Eq models with hq = 3. In all Eq models one expects 0J3 to run more slowly than in 
the MSSM due to the extra matter in the 5 and 5*' representations, keeping a 3 bigger and 
making two-loop contributions more important. We find that the latter destroy unification 
if there are three light generations of Eq matter (very similar to the case in Fig. |8|) even 
when an extra Higgs pair is included; there are model-dependent deviations from the curves 
in Fig. due to the presence of the extra U(l), but these are very small. Gauge coupling 
unification would require a^(Mz) to be reduced below the MSSM prediction by about 15%. 
Unification could conceivably still be rescued by threshold corrections from large splittings 
in the SU(5) multiplets at low energy. With three complete generations excluding the Hig- 
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gses, the situation is much worse as shown in Fig. |9|. 

(iv) SU(5) models with an extra 10 and 10*. This case is similar to case (iii), since the beta 
function coefficients 63 and 633 are the same. Two-loop contributions to the RGE's become 
relatively more important, and gauge coupling unification becomes problematic without 
large threshold corrections. 

We have followed a philosophy of preferring the least number of split representations 
possible. It is possible to make gauge coupling unification work without an intermediate 
scale far removed from the electroweak scale by relaxing this constraint. In fact, a non- 
super symmetric left-right Eq model has been proposed recently |56| in which the SU(2)r is 
broken at 1 TeV. 

(e) Yukawa evolution 



One can consider the evolution of Yukawa couplings in this new scenario where the QCD 
coupling does not run as rapidly as it does in the MSSM. For 63 = one can immediately 
solve the one- loop renormalization group equations (neglecting the SU(2) and U(l) couplings 
which are small except near the GUT scale). 

Consider the superpotential 



W = X t H 2 , 3 Qt c + \ b H h3 Qb c + X T H lt3 Lr c + X Si S 3 H 1:i H 24 + X di S 3 x di x c di 



(108) 



We define #1,3, #2,3 and S3 to be the linear combination of the Higgs doublets and singlets 
that acquire a vev. For the top Yukawa coupling one has (assuming that and A r can be 
neglected) 



dXt 



Xt 



Y / c l g 2 -3g 2 2 -^-g 2 3 + 6X 2 + X 2 S3 



(109) 



dt 16tt 2 

where ^ are the running U{1) gauge couplings (and hence contain some model dependence). 
The new couplings that arise from the presence of an electroweak Higgs singlet S evolve as 



dX Ss 


Xs 3 


dt 


16vr 2 


dX Sl 


x Si 


dt 


16tt 2 



J2d^-3g 2 2 + 3X 2 t +AX 2 S3 + 3j2Xl 



£ d l9 2 - 3g 2 2 + 3A 2 + 4X% 



i ^ 3 



(110) 
(111) 
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v, 



dt 167T 2 



i 6 j 



(112) 



where q, di and are some (in general model-dependent) coefficients of the U(l) gauge 
couplings. We do not assume any SU(5) relation between As^ and A^. 

Since the gauge coupling values near the GUT scale are much larger than they are in the 
minimal supersymmetric model, one expects the top Yukawa coupling to be driven much 
faster to its fixed point value from below. The general result is that if singlet quarks are 
accompanied by other exotics fermions to fill out representations of SU(5), then the top- 
quark is driven to its fixed point value (A^ ~ 8/9g%) over a wide range of values for the top 
quark GUT scale Yukawa; see Fig. [U]. The presence of the coupling X s in the top quark 
coupling RGE could soften the attraction to the fixed point, but requiring it to remain 
perturbative up to the GUT scale prevents it from destroying the fixed point solution, as 



shown in the NMSSM model 57 



The linear combinations H 13 , if 2 ,3 and S 3 of Higgs fields acquire vevs v 2 , t>i and v s and 
one defines tan/? = v-ijv\. One gets the usual relations that one has in the MSSM model 

m b (m b ) . . . V2m T (m T ) \/2m t {m t ) 

h(m t ) = w- , K(m t ) = a- , X t (m t ) = — — , 113 

r] b v cos p 7] T v cos p vsmp 

in addition to mass relations for the squark singlets and exotic leptons 

Ad <-^7' A ^"^7' (114) 

ignoring mixing. 



The singlet quark Yukawa couplings A^ also have infrared fixed points ||58|| , essentially 
given by the condition 

Unfortunately this does not yield a prediction for the singlet quark mass since the Higgs 
singlet vev is a priori unknown. At best one can obtain an upper limit on the ratio m Xd /Mz' 
5"8fl . When the singlet quark Yukawa is at its fixed point, it saturates this upper limit. 
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One can also consider the implications of bottom-tau unification in the context of these 
Eq models. The evolution of the Yukawa couplings is the same as it is for the MSSM, 



dR b / T 
dt 



R 



b/r 



167T 2 



4 2 16 2 i \2 , o \2 o\2 

-3^1 - y03 + x t +3A fe -3A r 



where Rb/ T = y-- However, we find that since the gauge couplings are larger over the entire 
range of scales between M GUT and the electroweak scale, the Yukawa couplings have to 
be correspondingly larger to cancel off the contributions from the gauge couplings. In the 
MSSM, the top Yukawa is often forced into the infrared fixed point region. In the E$ model 
with three light generations, we find that there is no solution that gives an acceptable value 
for m h f5"5"l- 



IX. CONCLUSIONS 

Quark singlets offer an interesting example of physics beyond the SM. They mix with 
the ordinary fermions. They impact a wide variety of experimental measurements, as they 
generate tree-level FCNC's, introduce unitarity violation in the SM CKM matrix, influence 
neutral meson-antimeson oscillations, and modify CP asymmetries. These objects can be 
produced by strong, electromagnetic and weak-neutral-current interactions, and produce 
interesting decay signatures. Their masses must generally exceed \Mz\ higher limits 85-131 
GeV apply in certain particular scenarios (see Section IV). 

We have collected the available bounds on singlet quark mixing; some have been up- 
dated; some, such as the B°, D° — > and D — > fi + [i~X constraints and the unitarity 
implications of FDNC bounds, have not appeared explicitly before (see Sections III- VI). 
The present limits on the 4x4 mixing matrix elements connecting one new singlet quark 
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to standard quarks may be summarized as follows. 
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We have discussed possible effects of singlet quarks on b — > e^y, 57 decays, and have 
illustrated how a w-type singlet could either increase or decrease the SM rate for b — > 37 
(Fig.|j). We have pointed out that small x — q mixing reduces the branching fraction for 
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Z — > qq decays; in the interesting case Z — > 66, this would worsen the present discrepancy 
between SM and experiment. We have given new illustrations of ways that singlet quarks 
can cause substantial deviations from SM expectations for CP-asymmetries of neutral B 
decays (Figs.||-|6]). The asymmetry Im X(Bd — > ipKs) can have the opposite sign to the SM 
value. 

In the GUT context, singlet quarks cannot simply be added by themselves to the SM or 
MSSM, since this would destroy gauge coupling unification (Fig.0); they must be accompa- 
nied by other members of exotic fermion multiplets. Down-type singlet quarks are readily 
accommodated in grand unified extensions of the SM; as a minimal scenario, they can be 
realized by adding one or more extra generations of 5 and 5* representations of SU(5), that 
imply extra vector-doublet leptons too. This exotic matter together with the SM matter 
content fits into the 27 representation of E 6 (which decomposes to 10 + 5* + l + 5 + 5* / -|-l in 
an SU(5) subgroup). Adding extra complete multiplets of SU(5) preserves (at the one- loop 
level) the successful unification of gauge couplings in the MSSM, since a complete multiplet 
contributes equally to the evolution of each coupling. However, more than three generations 
of exotic matter will destroy asymptotic freedom for a 3 at one loop (Fig||). 

As in the MSSM, gauge coupling unification in a desert model can be achieved by as- 
suming that split representations exist. In the context of models with singlet quarks, this 
means that there must be an additional pair of light Higgs doublets, in addition to the pairs 
that are included with the singlet quarks in the 5 and 5* representations. (The MSSM is 
then a special case consisting of no singlet quarks and one light pair of Higgs doublets). 

An up-type singlet quark is not contained as elegantly in GUT Models; it does not appear 
in the smallest representations, and its role is less clear. As a minimal prescription, it can be 
introduced by adding one extra light 10 and one 10* representation of SU(5) that get their 
mass from an SU(5) singlet Higgs boson; this implies extra vector-doublet quarks and a 
vector- singlet charged lepton too, preserving MSSM gauge coupling unification with b 3 = 
at one loop (Fig.^). Less minimally, it can also be realized in the SO (10) group with an 
extra light 45 (adjoint) representation (which decomposes to 24 + 10* + 10 + 1 in an SU(5) 
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subgroup), but this leads to nonperturbative gauge couplings at the GUT scale if the entire 
45 is required to be light. 

Two-loop effects are typically small in most GUT models, but if one includes extra 
representations of matter then the evolution of the strong coupling is diminished and it 
might even increase (no asymptotic freedom) toward the GUT scale. In a situation where 
the strong gauge coupling does not evolve at the one-loop level, we find that the two- 
loop effects become relatively more important and can make gauge coupling unification 
problematic, e.g. for three complete generations of E§ 27-plet matter. However one expects 
the low-energy threshold corrections to be more significant in this case, and ultimately the 
success of unification depends on the details of the low-energy spectrum. Two-loop effects 
also threaten perturbativity and asymptotic freedom of a 3 (Figs.|-|9|). 

Fixed points play a role in these extended model, with the top quark and the down-type 
singlet (s) masses possibly determined by the gauge couplings and the associated vevs. With 
extended matter content and larger gauge couplings, the top Yukawa coupling is driven 
to its fixed point faster than before (Fig.^). However, the Yukawa unification condition 
Xb{M G ) = X T (M G ) becomes harder to accomodate, and fails in the E 6 model with three light 
generations. 



X. APPENDIX A 

In this appendix we collect a few results needed for the analysis of b — > (q = s, d) 
in Section VI. The SM magnetic and chromomagnetic couplings for flavor-changing b —>■ q 
decays via W loops are given by [[59 



4'= 36(s- 1)' + iW' nI ' <116) 

= - ■ (ll7) 

where x = m^/M^. A down-type singlet quark induces additional Z loops, giving the 
replacements 
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These substitutions then enter into the values of the Wilson coefficients where 

"3 



c 7 (M w ) - 
c 8 (M w ) 



— -— sin 2 9 



V tb V t *J V54 81 



Zqb 



19 2 
4 ~ 81 
4 2 



2 Jg [ ' \V tb V t *J V9 27 



sin 2 # 



The coefficients from the 8x8 anomalous dimension matrix are IBS 



ai = ( i|, ±§, 0.4086, -0.4230, -0.8994, 0.1456) 

hi = ( §§f§, -§ff§, -f , -0.6494, -0.0380, -0.0186, -0.0057) 

In the case of up-type singlet quarks the Wilson coefficients are 



c\[M w ) 




c x 7 (M w ) 




c\{M w ) 


= \*f 


c x 8 (M w ) 


= lvfP{v), 



where y = m 2 Xu /M^. 



(120) 
(121) 



(122) 



(123) 
(124) 
(125) 
(126) 



XI. APPENDIX B 



We here collect some two-loop results needed in Section VIII. The two-loop RGE coeffi- 
cients for an arbitrary number of 10, 5*, and 5 representations are 



23 7 f N 9 

oil = jTj^io + (n 5 * + n 5 .# + n 5 ) + —n H 

3 9 9 

»12 = 77^10 + T77( n 5* + 7i5*' + ns) + , 
10 10 5 

6i3 = y^io + Y5^ 5 * + n5 *' + n ^ ' 



(127a) 
(127b) 
(127c) 
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13 3 
&21 = lO™ 10 + lO^ 5 * + n5 *' + n ^ + 5 Hh ' (127d) 
21 7 

^22 = y n io + ^ ( n s* + n 5*' + ^5) + 7n ff - 24 , (127e) 

& 23 = 8n 10 , (127f) 
3 2 

631 = 7^10 + Ys^ 5 * + n5 *' + n5 ) ' ( 127 g) 

632 = 3m , (127h) 

17 

633 = I7n w + —{n 5 * + n b *> + n 5 ) - 54 , (127i) 

in the g±, #2, g% basis. 

The two-loop coefficients for the E 6 models are 

1 1/ \ 81 25 97 

Kp = ^ n io + 7.K* + n 5 ) + —n 5 *, + —n N + —n H , (128a) 

3 1 9 39 

b ' P 2 = 2^ ni0 + s^ 5 * + n ^> + ^ n5 *' + Y00 UH ' ^ 128b ^ 

9 3 27 39 

frps = 2^ ni0 + 5^5* + n s) + ^q"™ 5 *' + ' ( 128c ) 

6 8 18 
frp 4 = 7^10 + -(n 5 * + n 5 ) + — n 5 *> , (128d) 

o 

„ 3 1. , 9 39 . . 

lp = 20 ni0 + 5^ 5 * + Ub > + 20™ 5 *' + T00 HH ' ^ > 

23 7 9 

6'n = 77^10 + ^( n 5* + n 5*' + n 5 ) + 7^n H , (128f) 

3 9 9 

K2 = 77^10 + 77(^5* + n 5 *> + n 5 ) + -ra H , (128g) 
1U 1U 

24 16 

&1 3 = -77^10 + 77 ( n 5* + ™5*' + n 5 ) , (128h) 
o id 

3 1 9 13 

b 2 P = 7^ n w + 5^* + n s) + 2Q™5« + , (128i) 

13 3 
& 2i = TTT^io + 77(^5* + n b *> + n 5 ) + -n H , (128j) 
1U 1U 

21 7 

b 22 = Y Uw + 2^ 5 * + U5 *' + 77,5 ^ _ 24 ' (128k) 

& 23 = 8nio , (1281) 

3 1 9 

6 ' 3p = 20 ni ° + o"^ 5 * + U ^ + 20 U5 *' ' (128m) 
3 2 

631 = 7^10 + 7^( n 5* + n 5*' + n 5 ) , (128n) 

6 32 = 3n 10 , (128o) 

17 

633 = 17ni + -77(^5* + ^5*' + ^5) - 54 , (128p) 
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10 



5 

16 4 . ,64 

y^io + 7~g( n 5* + n s*') + Y5 n5 ' 

2 1 8 17 
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o oU lo oU 
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Ki = TTT^io + ^( n 5* + n 5 *> + n 5 ) + — n H , (130f) 

1U oU ZD 

3 9 9 

b i2 = ttt^io + T7r( n 5* + n 5 */ + n 5 ) + -% , (130g) 

24 16 

b'" z = y^io + y^( n 5* + «5*' + n 5 ) , (130h) 

2 1 8 17 

&2p = g n io + ^7j( n s* + n 5") + + gQriif , (130i) 

13 3 
b 2i = TTT^io + T7:( n 5* + n 5*' + n 5 ) + -n H , (130j) 
1U 1U o 

21 7 

622 = y n io + 2 + "5" + "5) - 24 , (130k) 

b'£ = 8n 10 , (1301) 

2 1 8 

= ^10 + ^7j( n 5* + n 5*') + , (130m) 

3 2 

& 3 'i = -nio + Y5^ 5 * + n5 *' + n s) ' (130n) 

b'£ = 3n 10 , (130o) 
17 

633 = I7n 10 + — (n 5 » + n 5 */ + n 5 ) - 54 . (130p) 
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TABLES 

TABLE I. FCNC effects of singlet quarks 
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TABLE II. Quantum numbers of rank 5 Eq models 
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FIGURES 

FIG. 1. Unitarity of the 3x3 CKM matrix implies triangle relations like the one shown. In 
4x4 mixing cases the triangle relations become quadrangle relations. 

FIG. 2. The Z — > xx contribution to the toal Z decay width Tz is shown versus m x for 
Q x = -\ and §. 

FIG. 3. Singlet quark mixing can give meson-antimeson oscillations via induced FCNC tree 
diagrams (a) and via box diagrams (b), illustrated here for the B® — B^ d case. 

FIG. 4. Effects of a Q = | singlet quark on b — > decay. The branching fraction normalized 
to the SM value is shown versus the singlet quark mass. The curves are labelled by the values of 
|V^* S V^. 6 |; we assume that V* s V xb and V t *V tb have the same phase within an overall ± sign, shown 
on the label. 

FIG. 5. The CP asymmetry Im \{Bd — ► ipKs) in the presence of down-type singlet quarks. 
The band indicates the present uncertainty in the SM prediction. For some values of 5d and O^d 
(defined by Eq. fl83|)) there are no solutions as the unitarity quadrangle cannot be made to close. 

FIG. 6. The CP asymmetry Im \{B& 7r + 7r~) in the presence of down-type singlet quarks. 
The entire range is allowed in the SM. The parameters and 8^ are defined by Eq. ([83]). 

FIG. 7. One-loop gauge coupling evolution, adding either one d-type singlet (n Xd = 1) or one 
u-type singlet (n Xu = 1) to the MSSM. The singlets leave the SU(2)l gauge coupling unaffected 
at one-loop, but alter the evolution of the other gauge couplings and destroy unification. 

FIG. 8. One-loop and two-loop gauge coupling evolution with the addition of different numbers 
of light 10 multiplets of 50(10) to the MSSM. Successful gauge coupling unification is preserved 
with the addition of one or two 10-plets, but is threatened by two-loop effects when three light 
10-plets are added to the MSSM. 
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FIG. 9. Two-loop gauge coupling evolution where the electroweak scale matter consists of 
three light 27's. Gauge coupling unification is unsuccessful without the presence of a light pair of 
Higgs doublets from a split representation. The result for the case including a pair of light Higgs 
doublets is shown in Fig. (||). 

FIG. 10. The evolution of the top quark Yukawa coupling in the presence of three light 10 
multiplets of SO(10) added to the MSSM. The top quark Yukawa coupling reaches its infrared 
fixed point for a large range of initial (GUT) values, Xtc- The curves show \ t c between 0.2 and 4 
in increments of 0.2. 
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